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We present results showing that the strong coupling constant measured in two-flavor full QCD with dynamical 
Kogut-Susskind quarks at f3 = 5.7 exhibit a 15% increase due to sea quarks over that for quenched QCD at the 
scale fj, « 7GeV . 



1. Introduction 

Recently a lattice QCD determination of the 
strong coupling constant a a has been attempted 
employing the spin-averaged charmonium 15 — IP 
mass splitting to fix the scale[l]. The calculations 
are so far restricted to quenched QCD, and cor- 
rections due to sea quarks have to be estimated 
to make a physical prediction for a a . In this ar- 
ticle we report on our calculation of the strong 
coupling constant in two-flavor full QCD in the 
spirit of Ref. [1] and discuss implication of the 
results on the effects of sea quarks. 

2. Simulation 

Our study is based on the full QCD configura- 
tions on a 20 4 lattice previously generated with 
two flavors of dynamical Kogut-Susskind quarks 
at f3 = 5.7 with m q a = 0.01[2]. For charmo- 
nium spectrum measurement we employ a sub- 
set of these configurations, periodically doubled 
in the temporal direction. We use both Wil- 
son (without 'clover' improvement) and Kogut- 
Susskind actions for valence quarks. The Wilson 
results are obtained at K = 0.130 and 0.135 with 
72 and 75 configurations, and the preliminary 
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Table 1 

Charmonium operators. Smearing functions used 
for the Wilson case are f s {x) = exp( — |x| 2 /4) and 
P p (x) = sm(2irx>/ L)f s (x). e x = (-1)^+^3 

and Cxi = e x (-l) Xi . 

meson source sink 

(a) Wilson quark 

MP1) E xy ^<T ik rl>yfs(x)fs(y) $*<T> k rl> x 

MP2) J2 x , v i>x-Y 5 i> v fU x )fs(y) f^jh 

(b) Kogut-Susskind quark 

J /i> Y<xCxiXx,tXx,t 

Vc(NG) J2x e *Xx,tXx,t 

VC J2 X (Xx,tXx,t + l - X.X,t + lXx,t) 

Xx,t+1 + Xx,t+lXx,t) 



Kogut-Susskind results, obtained after the confer- 
ence, are based on 19 configurations at m q a = 0.3 
and 0.4. 

For Wilson valence quarks we used a Gaussian 
smeared source and a local sink. Typical forms 
of the operators are listed in Table 1. Two types 
of operators are used for the 1 Pi state h c , labeled 
PI and P2. 

For Kogut-Susskind valence quarks we employ 
operators constructed in Ref. [3], which are also 



Table 2 

Charmonium 15 



IP splitting in full QCD. 



operator K/m q a 


Am 1S -ipa 


7r/a(GeV) 


(a) Wilson valence 


quark 




/i c (PI) 0.130 


0.221(49) 


6.5(1.5) 


0.135 


0.221(48) 


6.5(1.4) 


h c (P2) 0.130 


0.200(32) 


7.2(1.2) 


0.135 


0.205(35) 


7.0(1.2) 



(b) Kogut-Susskind valence quark 

r] c {NG) 0.3 0.185(20) 7 

0.4 0.208(23) 6. 

rj c 0.3 0.148(20) 9. 

0.4 0.164(23) 8. 



79(84) 
93(77) 
7(1.3) 
8(1.2) 



listed in Table 1. For rj c the first operator corre- 
sponds to the Nambu-Goldstone channel of U(l) 
chiral symmetry of the Kogut-Susskind action. 
For sources we employ an improved wall source 
technique described in Ref. [4]. 

3. Determination of scale 

3.1. Wilson valence quarks 

In Table 2(a) we list our result for the charmo- 
nium 15 — IP mass splitting obtained with Wilson 
valence quarks and the corresponding scale ir/a. 
We observe unusually large errors of 20-25% in 
spite of the use of over 70 configurations. In con- 
trast a test run in quenched QCD on a 16 3 x 32 
lattice at f3 = 6.0 with 20 configurations showed 
only a 7% error, and the full QCD results with 
Kogut-Susskind valence quarks described below 
also exhibit smaller errors with 19 configurations. 
A possible cause of the large errors for Wilson va- 
lence quarks might be ascribed to the mismatch 
of quark actions taken for sea and valence quarks, 
though actual connection is not clear to us. 

3.2. Kogut-Susskind valence quarks 

It is more natural to employ the Kogut- 
Susskind quark action for valence quarks on our 
full QCD configurations generated with the same 
action for dynamical sea quarks. Our results for 
the 15 — IP splitting obtained for this case are 
listed in Table 2(b). 

In Fig. 1 we show the charmonium spectrum 
for the valence quark mass m q a = 0.3. Con- 



> 
a) 

a 35 



KS valence quark ( m a=0.3 1/a=2.48(27)GeV ) 
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Figure 1. Charmonium spectrum in two-flavor 
full QCD at (3 = 5.7 with Kogut-Susskind sea 
(m q a = 0.01) and valence (m q a = 0.3) quarks. 
Errors shown do not include that of scale. Hori- 
zontal lines indicate experimental values. 



version to physical units is made with ir/a ob- 
tained with 7j c (NG) given in Table 2. Except 
for the degeneracy of J/ijj and 7j c in the non- 
Nambu-Goldstone channel, the pattern of the 
spectrum is consistent with the experiment. In 
particular 7j c in the Nambu-Goldstone channel is 
lighter than J ftp (the value of the mass splitting 
m J/ij> ~ ™rj c = 360(40)MeV is three times larger 
than the experimental value 118MeV, however). 
This is quite different from the results with the 
Wilson action for which the addition of the clover 
term is needed to split the two 15 states. We note 
that the errors in the physical values of masses 
shown in Fig. 1 do not include that of the scale 
7r / a of 11%, which translates to about 400MeV 
for the masses. An increase in statistics is needed 
to see if the physical charm quark corresponds to 
a smaller value of m q a. 

4. Strong coupling constant 

Our full QCD results for the 15 — IP mass split- 
ting yield a value ir/a « 7GeV with an error of 
10-20% for Wilson and Kogut-Susskind valence 
quarks (the latter with the Nambu-Goldstone 7j c ). 
This value is consistent with ir/a = 7.01(28)[2] 
estimated from the p meson mass using Kogut- 
Susskind valence quarks. The use of the non- 
Nambu-Goldstone operator for 7j c , however, leads 
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Figure 2. Evolution of a g (p) for quenched and 
two-flavor QCD. Thick bars show measured val- 
ues of a g (p) and uncertainties in the scale. 



to a value 30% larger for the scale. This discrep- 
ancy represents a systematic uncertainty for the 
Kogut-Susskind case. 

In order to estimate the strong coupling con- 
stant at the scale estimated above, we employ the 
relation[l,5], 

ai^Tr/a) -1 = Pa _1 + 0.3093 - 0.0885JV/, (1) 

where the last term represents the contribution of 
Nf flavors of Kogut-Susskind sea quarks[6]. With 
P = 0.577323(34)[2] at /3 = 5.7 and N f = 2 this 
gives 

c£L = 0.142. (2) 

MS v > 

For comparison we quote the quenched result: 
a£lL = 0.124, (3) 

MS ' v ' 

for ir/a = 7.63(48)GeV at (3 = 6.1[1]. We thus 
find that the strong coupling constant extracted 
for two-flavor full QCD is about 15% larger than 
in the quenched case at fi « 7GeV. 

A larger value is in fact expected from the con- 
sideration that fixing the scale by the IS — IP 
mass splitting means adjusting the QCD cou- 
pling strength at a scale typical of charmonium 
and that the coupling for larger momenta de- 
creases more slowly in the presence of sea quarks. 
This view underlies the heuristic procedure of 



Ref. fll in which the value of a—— was esti- 

L 1 MS 

mated from the measured value of ai— L by match- 

MS J 

ing the two couplings evolved down to the scale 
H « 0.4 — 0.75GeV with the two-loop renormal- 
ization group (3 function. 

(2) 

In Fig. 2 we apply this procedure to otjjg and 
ai^L starting with (2) (the Wilson result with 
h c (P2) at K = 0.135 is used for 7r/a; other cases 
are similar) and (3). The overlap of the two 
bands at fi « 0.5GeV shows that a 15% increase 

(2) 

found for a—— is quantitatively consistent with 
the expectecfmagnitude of sea quark effects. This 
agreement also implies that our two-flavor result 
would lead to a value of ot—^— similar to the origi- 

MS ° 

nal result allL.(5GeV) = 0.172(12)[1]. In fact we 
find allL.(5GeV) = 0.172±°;°^ usmg our Wilson 
result with h c (P2) at K = 0.135 in Table 2. 

5. Summary 

Our results and analyses show that sea quark 
effects are visible in the strong coupling constant 
measured in current full QCD simulations. This 
indicates a promising prospect in the near future 
for a full QCD determination of the strong cou- 
pling constant including a more realistic spectrum 
of sea quarks than was attempted here. 
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